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In helical hydromagnetic turbulence with an imposed magnetic field (which is constant in space 
and time) the magnetic helicity of the field within a periodic domain is no longer an invariant of 
the ideal equations. Alternatively, there is a generalized magnetic helicity that is an invariant of the 
ideal equations. It is shown that this quantity is not gauge invariant and that it can therefore not 
be used in practice. Instead, the evolution equation of the magnetic helicity of the field describing 
the deviation from the imposed field is shown to be a useful tool. It is demonstrated that this tool 
can determine steady state quenching of the alpha-effect. A simple three-scale model is derived to 
describe the evolution of the magnetic helicity and to predict its sign as a function of the imposed 
field strength. The results of the model agree favorably with simulations. 

PACS numbers: 52.65.Kj, 47.11. +j, 47.27.Ak, 47.65.+a 
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I. INTRODUCTION 

Magnetic helicity has traditionally been used as a diag- 
nostic tool to characterize magnetic field topology. Only 
in recent years magnetic helicity has also become a useful 
tool in understanding large scale dynamo action. Mag- 
netic helicity is important because it is conserved in the 
limit of vanishing resistivity. This is not the case with 
the kinetic helicity, which is also conserved in the inviscid 
case, but the kinetic helicity dissipation rate diverges in 
the inviscid limit 0. In this sense kinetic helicity is not 
even approximately conserved at large Reynolds num- 
bers, while magnetic helicity is very nearly conserved at 
large magnetic Reynolds numbers @]. 

A key result that has emerged from the concept of 
magnetic helicity conservation is that, in a periodic do- 
main, a large scale magnetic field generated by the a- 
effect 0, saturates on a resistive time scale [|| . This 
time scale can be very long. The helicity concept has also 
provided us with a simple explanation for the final sat- 
uration field strength of helical dynamos in periodic do- 
mains; see Ref. [j| for details. In this light, the case with 
an imposed magnetic field has also been considered 
where it was found that above a certain field strength the 
dynamo is suppressed. 

In the present paper we use similar ideas to obtain a 
more detailed understanding of the case with an imposed 
field. We begin with the equation governing the evolu- 
tion of the magnetic helicity in a periodic domain in the 
presence of an imposed field. It is well known that in 
that case the magnetic helicity of the fluctuating mag- 
netic field is no longer conserved in the nonresistive limit 
but we also point out that a certain generalized (or 
total) magnetic helicity that has sometimes been used in- 
stead is gauge dependent and can therefore not be used 
in the present case. We then discuss applications of the 
magnetic helicity equation to the alpha-effect in mean- 
field electrodynamics and derive a model equation in or- 



der to understand the evolution of the magnetic helicity 
for different imposed field strengths. 



II. MAGNETIC HELICITY EQUATION 

The evolution of the magnetic field B is governed by 
-VxE, V-B = 0, (1) 
where the electric field is obtained from Ohm's law, 



dB 

~dt 



E 



u x B + r)J, 



(2) 



where u is the velocity, J the current density, and 7] the 
resistivity. Throughout this paper we adopt SI units, 
but we set the permeability to unity. We consider all 
quantities to be triply-periodic over a cartesian domain. 
We consider the case with a finite mean field, 



(B) = B = const ^ 0, 



(3) 



where angular brackets denote full volume averages. Such 
averages have no spatial dependence, but they can still 
depend on time. However, because of periodicity, the 
volume average of the curl in Eq. vanishes, and hence 
d(B)/dt = 0. In other words, Bq is not only constant in 
space, but it is also constant in time. 

Next, we split the field into a mean and a fluctuat- 
ing component, B = Bq + 6, and introduce the mag- 
netic vector potential for the fluctuating component via 
b — V x a, where a is periodic. The uncurled induction 
equation reads 



da , _ 
where <f> is the scalar potential. 



(4) 
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We now consider the magnetic helicity. The imposed 
field is constant in space and does therefore not con- 
tribute to the magnetic helicity. We therefore con- 
sider only the magnetic helicity of the fluctuating field, 
H = (a ■ b). This quantity is gauge invariant because 
adding a gradient term to a does not change H: 

((a + V<p)-b)=H+{V- (ipb)) - {ipV b)=H. (5) 

Here we have used the solenoidality of b, and the fact 
that the volume average over a divergence term vanishes 
for a periodic domain. The equation for the (gauge- 
dependent) helicity density of the fluctuating field can 
be obtained in the form 

^(a-b) = -2£-6-V-p-V(/i)xo]. (6) 

The divergence term vanishes after volume averaging, so 

±{a-b) = -2{E-b), (7) 

where all terms are gauge-independent. Making use of 
Eq. J5J, we have 

±{a-b} = 2((uxB )-b)-2r 1 (j-b), (8) 

where we have used (J-b) = (j ■ b). Equation JSJ) can 
also be written as 

^-(a-b) = -2£ -B -2 V (j-b), (9) 
at 

where the electromotive force, £o — £o(t) = (u x b), has 
been introduced. If the flow is isotropic and helical, there 
will be an a-effect with £ = aB , so 

^ t (a-b) = -2aBl-2 V (j-b). (10) 

(In Appendix 1X1 we clarify the implications of a finite a 
effect when the Faraday displacement current is restored 
in the Maxwell equations and when it is ignored.) Since 
(a • b) is gauge invariant, it is a physically meaningful 
quantity. If there is a steady state, then (a • b) must also 
be steady. In that case we have 

a = -V(j ■ b)/B 2 , (11) 

which is a relation due to Keinigs [fj for the a-effect in 
the saturated (steady) state; see also Ref. ^(j- If t ne 
field is weak, the a-effect will remain finite in the high 
conductivity limit 0. 

The presence of a finite a-effect means that the struc- 
ture of Eq. I|ll|) is very different when there is an imposed 
field. Unlike the case without imposed field (Bo = 0), 
the quantity H = (a ■ b) is no longer conserved in the 
limit 77 — ^ 0. This prompted Matthaeus & Goldstein Q 
to consider the quantity 

H = H + 2A -B , (12) 



where 

A (t)= f £ (t')dt'. (13) 
Jo 

Note that H is constructed such that it satisfies the equa- 
tion 

— = -2 V {j-b). (14) 

Indeed, this equation reduces to Eq. @ after inserting 
Eqs (T3J| and into (Tg) . 

The symbol A is chosen to make the extra term in 
Eq. (|12fl look like a magnetic helicity, even though V x 
j4 = ^ B . The reason for V x A ^ B is that 
Aq corresponds to a slowly varying variable, whose curl 
gives Bq at a higher order, and not the order we are 
working in; see Eq. (A9) of Ref. 11]. Conversely, B 
corresponds to the curl of Aq at a lower order. A rigorous 
scale expansion is given in the Appendix of Ref. [IlJ • 

In the nonresistive limit, the right hand side of Eq. 1141) 
vanishes, and so H is conserved. One may be tempted 
to conclude that in the steady state, (j ■ b) = 0. This is 
not generally true, however, and it would be in conflict 
with Eqs l|10fl and (|llfl . Certainly for sufficiently weak 
fields a is finite 12], so rj(j ■ b) will also remain finite, see 
Eq. HI ip . Therefore, H cannot be constant in the steady 
state. The reason for this puzzle is that H is not gauge 
invariant 0], because the definition of H involves the 
quantity A . At first glance, A appears to be gauge 
invariant, because £q is gauge invariant and Aq involves 
only a time integral over £q; see Eq. (|13|) . However, 
the beginning of the time integration is ill-defined, so in 
general one can replace 

A (t) -» A Q (t) = A 00 + f £ (t') dt', (15) 
Jo 

which would lead to a different conserved magnetic he- 
licity, H + AH, where AH = 2Aoq ■ Bq = const is un- 
determined. Therefore, H is not a physically meaningful 
quantity, so it is not surprising that H can have a compo- 
nent that grows linearly in time. We emphasize however 
that H = (a ■ b) is still gauge invariant and therefore 
physically meaningful, even though it is no longer a con- 
served quantity in the ideal limit. 

III. NON-PERIODIC GAUGE POTENTIALS 

In this section we want to comment on the related is- 
sue that adding a spatially constant vector Eo(t) to the 
right hand side of Eq. (0J would not affect the evolution 
of b. The constant vector E {t) can readily be absorbed 
in the definition of the scalar potential </>, because it is 
specified only up to an additional gauge potential. How- 
ever, this gauge potential has in general a non-periodic 
contribution, even when all other quantities are periodic. 
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More specifically, (j) in Eq. @ must have an additional 
component that varies linearly in space, i.e. 

4> = 4>-E -x, (16) 

where <j) is periodic and x is the position vector. We stress 
that 4> in Eq. 1)16(1 is therefore in general not periodic, 
even if a and b are periodic. 

A comment on the helicity flux associated with the 
gauge field is here in order. This flux is often written 
as 4>b which would then not be periodic and hence it is 
not obvious that its surface integral vanishes; see Eqs © 
and Q) . (The importance of this term for magnetic he- 
licity injection has been discussed in Ref. However, 
using the identity 

V<f) x a = V x (cj)a) - cf)b, (17) 

the magnetic helicity flux 4>b can also be written as V</> x 
a, which is periodic. Therefore, there is no contribution 
from the (j)b term in our case. (We note that similar 
manipulations can be used to turn a non-periodic vector 
potential into a periodic one if the velocity is a linear 
function of coordinates 0|.) The term V x (4>a) has 
recently been discussed in a formulation of a magnetic 
helicity conserving dynamo effect [16j . Obviously, such a 
term does not give a contribution under the divergence 
and hence cannot be physically meaningful 



IV. APPLICATION TO a-QUENCHING 

There have been a number of simulations of helically 
forced periodic flows with an imposed magnetic field. 
The general objective is to obtain the a-effect and its 
suppression as a function of field strength 0, 0, 0] . 

In the steady state, Eq. (fTT|l can be used to determine 
a by measuring (j ■ b) in a simulation with an applied 
magnetic field Bo- For helical turbulence, (j ■ b) can be 
approximated by 

(j-b) «e f fc f (b 2 ), (18) 

where et — ±1 for a fully helical field with positive or 
negative helicity, and |ef| < 1 for fractional helicity. A 
strongly helical small scale magnetic field is generally ex- 
pected when the turbulent velocity field is also strongly 
helical [2(J . The steady state a is therefore given by 

a = — €fkff](b 2 )/Bl (steady state value). (19) 

On the other hand, the kinematic value of a can be esti- 
mated in terms of the kinetic helicity, 

oik — — |t(o; • u) (kinematic value), (20) 

where r is the correlation time which, in turn, can be 
expressed in terms of the turbulent magnetic diffusivity 



for which we have a similar expression, r} t — jt(u 2 ). In 
analogy to Eq. I|18|l . we write (u> • u) » effcf (u 2 ), so 

a fi (b 2 ) 

— w — 7t (steady state value). (21) 

a K Vt ' ' 

If we assume that the small scale field is in equipartition, 
i.e. (b 2 ) « (/dopu 2 ) = B 2 q , and if we define [2l| the mag- 
netic Reynolds number as R m = rj t /rj, then Eq. I|21|) can 
be turned into the interpolation formula, 

a= l + R TB 2 /B^ (22) 

that recovers Eq. 1(21(1 for strong fields and a = q;k in 
the weak field limit, Bq — * 0. This equation is known 
as the catastrophic quenching formula of Vainshtein and 
Cattaneo (l2|. 

In order to confirm that the onset of steady state 
quenching depends on the magnetic Reynolds number 
based on the forcing scale j^, R m ,f = u rms /(?]kf), and 
not on the magnetic Reynolds number based on the scale 
of the box |3, R m ,i = u ms / (rjki), we show two series 
of simulations obtained for different values of the forcing 
wavenumber kf, for different values of Bo. The forcing 
of the flow was fully helical; for details on the numerical 
method see Ref. 5] . The result is shown in Fig. ^ 

Next, we consider simulations where we use hyperdif- 
fusivity, i.e. the ordinary magnetic diffusion operator, 
r]V 2 B, is replaced by (— l) n ~ 1 r] n V 2n B, where n = 1 
corresponds to the standard case. This is a common tool 
in order to extend the inertial range of the turbulence 
[23|. but it is also clear that this leads to wrong sat- 
uration field strengths J24J. In the presence of hyper- 
diffusivity, the magnetic Reynolds number is defined as 
/(rink 271 ' 1 ). In the following we show, how- 
ever, that the quenching data are better described by a 
single quenching curve when fcf is rescaled, 

-R m ,fs = Urms/^nfcf" -1 ), where fc fs = 1.6 X kf. (23) 

The result is shown in Fig. [5] The factor of 1.6 is not 
expected to be universal but is probably a slowly varying 
function of magnetic Reynolds number |2lL l24j . 

We conclude that Eq. l|22^l describes the simulations 
quite well provided the magnetic Reynolds number is de- 
fined in a suitable manner. We emphasize however that 
this equation only applies to the steady state and if there 
is no mean current. This is generally not the case and 
therefore the quenching is in practice not automatically 
catastrophic |2j], i.e. the onset of quenching docs not 
depend on R m . 

V. EVOLUTION OF LARGE SCALE 
MAGNETIC HELICITY 

Recently, the effect of an imposed field on the inverse 
cascade has been studied 0, 0- ^ tne imposed mag- 
netic field is weak or absent, there is a strong nonlocal 
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FIG. 1: Normalized a-effect versus normalized magnetic en- 
ergy, scaled with the small scale magnetic Reynolds number 
(upper panel) and the large scale magnetic Reynolds number 
(lower panel). Note that the onset of quenching is governed 
by the small scale magnetic Reynolds number, not the large 
scale magnetic Reynolds number. 



transfer of magnetic helicity and magnetic energy from 
the forcing scale to larger scales. This leads eventually 
to the accumulation of magnetic energy at the scale of 
the box 0, 123L |25| . As the strength of the imposed field 
(wavenumber k — 0) is increased, the accumulation of 
magnetic energy at the scale of the box (k = 1) becomes 
more and more suppressed 0. 

Qualitatively, this can be understood as the result of 
two competing effects: (i) the inverse cascade that pro- 
duces magnetic helicity of opposite sign at k = 1 com- 
pared to that at the forcing wavenumber kf , and (ii) the 
a-effect operating on the imposed field producing mag- 
netic helicity of the same sign at k = 1 than at k = kf. 
This is because the sign of the a-effect is opposite to the 
sign of the magnetic helicity at k = fcf, and a enters 
with a minus sign in the evolution equation (|10|) of mag- 
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FIG. 2: Normalized a-effect versus magnetic held strength. 
The full dots denote runs where hyperdiffusion has been used. 
In the second panel, Bq/B^ is scaled with the magnetic 
Reynolds number based on the forcing scale, Rm,f- The 
last panel is similar to the second, but the effective forcing 
wavenumber has been used which is kf scaled by a factor 1.6. 
This brings especially the hyperdiffusive runs (full dots) closer 
to the rest of the data points. 



netic helicity. Under the assumption that the turbulence 
is fully helical, the critical value of the imposed field 
can be estimated by balancing the two terms on the right 
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FIG. 3: Evolution of the total magnetic helicity, H = Hx+Hf, 
as a function of t for different values of Bo, as obtained from 
the three-dimensional simulation. 



hand side of Eq. i|l(J|) and by approximating, as in 5 IIVI 
a ~ Vt e fk{ and (j ■ b) ps efkfB^ . This yields 



£?/B e 2 q « V/Vt = R n 



(24) 



where the last equality is again to be understood as a 
definition of the magnetic Reynolds number, see also 
Ref. [2]]]. For Bq > the sign of the magnetic he- 
licity is the same both at k — 1 and at k = kf , while for 
B < B v the signs are opposite. 

A related phenomenological model for saturation of the 
dynamo effect under influence of Bq has been given 
that is based upon a Fourier scale separation approach. 
That approach leads the conclusion that the critical value 
£?* oc i?" 1 rather than Rm 1 ^ 2 as above. Further analysis 
may be needed to fully reconcile the differences in these 
approaches, both of which appear to have some support 
from simulations. 

A more quantitative description of the evolution of the 
magnetic helicity can be obtained by using a modified 
two-scale model [21], |23 , where the term 2£q ■ Bq from 
Eq. JjJJ has been included, so 

Hi = -2rjk 2 Hi + 2{£ 1 ■ B x ) - 2£ ■ B , (25) 



H i = -2rik%Hf-2{e l -Bi}. (26) 

Here, Hi and Hf are the magnetic helicities at the 
wavenumbers 1 and kf, respectively, and (£1 • B\) is the 
helicity production from a-effect and turbulent diffusion 
operating on the field at k = 1. We note that the sum 
of Eqs and yields Eq. Q. The electromotive 
force £1 at wavenumber ki is given by 



£1 = a,{Bx - rj t Ji 



(27) 



To calculate {SfBi) in Eqs and ® we dot Eq. (57)) 
with £?i, volume average, and note that (Ji-Bi) = k\Hi 



and (-Bf) = The latter relation assumes that the 

field at wavenumber ki is fully helical, but that it can 
have either sign. Thus, we have 

(£ 1 -Bi)=a { ki\H 1 \-r l tkiHi. (28) 

The large scale magnetic helicity production from the a- 
effect operating on the imposed field is £q ■ Bq = uiB\. 
The a-effect is proportional to the residual magnetic he- 
licity of Pouquet, Frisch and Leorat |2£j, with 



a = -\t{{ui ■ u) - (j ■ b)/p ) , 



(29) 



where r is the correlation time and po the average density. 
In terms of Hi and Hf we write 



*i = a K + hrklH 1} 



Otf 



= a K + WkfHf, 



(30) 



(31) 



for the a-effect with feedback from Hi and Hf, respec- 
tively. Here, qk is the contribution to the a-effect from 
the kinematic helicity, as defined in Eq. (|20|l . 

The above set of equations for the case of an imposed 
magnetic field is similar to a recently proposed four-scale 
model 26], where two smaller scales were added relative 
to the two-scale model. In the present case, on the other 
hand, instead of including scales smaller than the forcing 
scale, the imposed field at the infinite scale is included, 
albeit fixed in time. 

For finite values of Bq, the final value of Hi is partic- 
ularly sensitive to the value of ax and turns out to be 
too large compared with the simulations. This disagree- 
ment with simulations is readily removed by taking into 
account that q;k = — |r(w • u) should itself be quenched 

Thus, we write 



when Bq becomes comparable to B cq 
Q!ko/(1 



a K 



B o/ B cq)i 



(32) 



which is a good approximation to more elaborate expres- 
sions [53 ■ We emphasize that this equation only applies 
to Qk and is therefore distinct from Eq. (|22|l . 

In Fig. 01 we show the result of a numerical integra- 
tion of Eqs (|25|l and l|26|) . Both the three-dimensional 
simulation and the two-scale model show a similar value 
of Bq « 0.06. ..0.07, above which Hi changes sign. This 
confirms the validity of our estimate of the critical value 
£?* obtained from Eq. I|24|) . Secondly, the time evolution 
is slow when Bq < B* and faster when Bq > B*. In the 
simulation, however, the field attains its final level for 
Bq > B* almost instantaneously, which is not the case in 
the model. It is possible [7j that the almost instantaneous 
adjustment in the simulations is a consequence of the 
Alfven effect, which is not included in the present model. 
This, and other shortcomings of the present model may 
also be responsible for the mismatch between the mag- 
netic helicity amplitudes seen in the simulations and the 
model. Most characteristic in the simulations is the fact 
that Hi — > while Hf ^ in the limit of strong imposed 
field strength. 
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FIG. 4: Evolution of magnetic helicity as a function of t for 
different values of Bo, as obtained from the two-scale model. 



VI. CONCLUSIONS 



|26|. The nonlinear two-scale and multi-scale models 
play important roles in modern mean-field dynamo the- 
ory. Given that we are still lacking a proper understand- 
ing of solar and stellar dynamos in the nonlinear regime, 
an independent confirmation of the nonlinear multi-scale 
model must therefore be regarded as a crucial step toward 
understanding the origin and maintenance of magnetic 
fields in turbulent astrophysical bodies. 
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We have shown that (i) in the presence of an imposed 
field and (ii) using triple-periodic boundary conditions, 
the generalized magnetic helicity Js] in Eq. I|12fl is not 
gauge-invariant and can therefore not be used for prac- 
tical purposes. This quantity has frequently been used 
in the solar wind community as an alternative to the 
ordinary magnetic helicity which is known not to be con- 
served in the limit of vanishing resistivity. We have ar- 
gued, however, that even through the ordinary magnetic 
helicity is not conserved in the presence of an imposed 
field and in the limit of low resistivity, it remains an ex- 
tremely useful quantity that has predictive power - sim- 
ilar to the case without imposed field Q. 

Based on analytic considerations and confirmed by the 
simulations, we have shown that the sign of the magnetic 
helicity depends on the strength of the imposed magnetic 
field. If the field is weak enough, the situation is similar 
to the case without imposed magnetic field and the sign 
of the magnetic helicity is opposite to the sign of the 
helicity of the turbulence. If the field exceeds a certain 

— 1/2 

threshold, which is R m times the equipartition field 
strength, where R m is the magnetic Reynolds number 
based on the forcing wavenumber, the sign of magnetic 
helicity changes and becomes equal to the sign of the 
helicity of the turbulence. This can be understood as 
a consequence of the a-effect operating on the imposed 
field. In finite systems, this a-effect would cause the large 
scale field to have opposite helicity compared to the small 
scale field. In an infinite (or periodic) system, this is not 
possible, and the entire field in the computational domain 
plays the role of a small scale field which must then have 
the same sign of helicity as the turbulence. 

The two-scale model used to describe the nonlinear 
evolution of helical dynamos 0, can be generalized 
to take account of the large scale field. The formal- 
ism is similar to a recently proposed four-scale model 



APPENDIX A: THE ROLE OF THE 
DISPLACEMENT CURRENT IN TRIPLY 
PERIODIC SYSTEMS WITH IMPOSED FIELD 

In this appendix we discuss a paradoxical situation 
|28j that arises when comparing volume averages of the 
present equations (where periodic boundary conditions 
are used and a uniform field is imposed) with the full 
Maxwell equations (where the displacement current is in- 
cluded). The displacement current is given by the fourth 
Maxwell equation, 



1 dE 



= V x B J. 



(Al) 



where c is the speed of light. We recall that the perme- 
ability has been put to unity. Applying volume averages, 
and noting that the volume average of the curl vanishes, 
we have 



c 2 dt ~ { h 



(A2) 



On the other hand, the volume average of Ohm's law (J2J 
yields 



(E) = -(« x B)+r](J). 



(A3) 



We assume that there is no net flow, i.e. (u) = 0. There- 
fore, (u x B) = (u x 6). 

In helical hydromagnetic turbulence there is an a-effect 
H H so S = (u x b) = aB ^ 0, and therefore 
(E) ^ and, because of Eq. JUJl, (J) ^ 0. The latter- 
condition is, of course, inconsistent with our assumption 
that (J) = (V x B) = 0. This discrepancy could be par- 
ticularly important when considering the contribution of 
the volume averaged field to the Lorentz force, (J) x (B), 
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TABLE I: Comparison of various volume averages in the pre- 
Maxwellian approximation and the full Maxwell equations. 
The asterisk denotes that this result is strictly valid only when 
also a time average is considered. 





pre-Maxwell 


Maxwell 


(u x b) 






(E) 






(J) 


= 




(J) x (B) 


= 


= 0* 



which vanishes in the pre-Maxwellian MHD approxima- 
tion, but not when the displacement current is retained 

m. 

Eliminating (E) from Eqs ljA2() and HA3fl we have 

(c 2 +^(J)=5 , (A4) 

where the dot on £q denotes time differentiation. Equa- 
tion (|A4|) could be solved for (J) either in terms of the 
Green's function exp[— (i — t')c 2 /rf\ or via series expan- 
sion [2^, confirming that («7) ^ 0. However, here we are 
interested in the Lorentz force, so we write 

^ 2 + ^j ({J)x(B))=f xBo = 0. (A5) 

The right hand side of Eq. I|A5() vanishes, because (B) = 
Bo is independent of time and £q — c\Bq is parallel to 
B . The use of the equation £ = aB ignores random 
fluctuations in time about zero. In that sense, Eq. (|A5|) is 
strictly valid only when the averages are also taken over 
time. 

The solution to Eq. l|A5j) shows that, if the Lorentz 
force from the mean field was vanishing initially, it must 
vanish at all times. Table [I] summarizes which of the 
different volume averages discussed in this section vanish 
in the pre-Maxwellian MHD approximation and which 



quantities remain finite when the full Maxwell equations 
are used. 

The apparent inconsistency is removed by noting that 
Eq. (|AI|) does simply not exist in the pre-Maxwellian 
MHD formulation and hence cannot be invoked in the 
discussion. [The situation is similar to the incompress- 
ibility assumption, V-tt = 0, or the anelastic approxima- 
tion, V- (pu) — 0, both of which do not imply dp/dt = 0. 
Indeed, the original continuity equation is no longer used 
and has instead been replacedhy Vm = or V-(ptt) = 0, 
respectively.] Nevertheless, as far as the Lorentz force is 
concerned, the neglect of the displacement current is in- 
consequential, because it vanishes in either of the two 
cases; see Table Q] 

The mismatch between (J) = in the pre-Maxwellian 
approximation and the exact result, (J) ^ 0, is negligi- 
ble, but can be quantified using a rigorous ex pan sion in 
terms of slowly and rapidly varying variables . Such 
an approach also demonstrates quite nicely that the diffi- 
culties introduced into the periodic model by the presence 
of a nonzero uniform mean field are due to imposing peri- 
odic boundary conditions on the entire (infinite volume) 
system. If instead (see Appendix of Ref. II]) the turbu- 
lence is assumed to be modeled as locally homogeneous, 
in the statistical sense, and periodicity is employed in a 
two scale expansion as a local leading order model, no 
such problems emerge. 

Paradoxical situations arising from the assumption of 
triple periodicity are commonly resolved using scale ex- 
pansion. Another such example is the famous Jeans swin- 
dle [29J, where the assumed zero order equilibrium state 
does not obey triple periodicity; see Refs [30, H2 for a sta- 
bility analysis using a proper equilibrium solution. We 
emphasize, however, that the problem with the Jeans 
swindle is distinct from the problem with the displace- 
ment current discussed here. The latter is completely 
resolved by staying fully within the pre-Maxwellian for- 
mulation, while the former is a true mathematical swin- 
dle. 
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